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p Abstract 



An improved volume-weighted probability measure for eternal inflation is proposed. 
For the models studied in this paper it leads to simple and intuitively expected gauge- 
invariant results. 



O ■ 1 Introduction 

o 



Soon after the invention of inflationary cosmology it was realized that inflation may divide our 
universe into many exponentially large domains corresponding to different metastable vacuum 
states, forming a huge inflationary multiverse [TJ [2] . The total number of such vacuum states 
in string theory can be enormously large [21 IH El E]- A combination of these two facts with the 
eternal inflation scenario [Tj, [TJ [S] is called the theory of the inflationary multiverse [HI QUI EH], 
or the string landscape scenario [T2] . 

In order to make any kind of predictions in the context of this new scientific paradigm, 
one should find all possible vacua of the theory, study all possible versions of the cosmological 
evolution which can bring us to each of these vacua, and assign certain probabilities to all 
possible outcomes of this evolution compatible with the fact of our existence. 

This last part of the problem is especially complicated and ambiguous; see [131 [EU d5] for 
a recent discussion. For example, one may try to calculate the probability to be born in a 
part of the universe with given properties at a given point, ignoring creation of new volume 
(new points) during inflation. One can consider a small box of a given size Ax in comoving 
coordinates, which do not change during the expansion of the universe, take a representative 
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set of points in this box, and follow the evolution of physical conditions (fields, temperature, 
etc.) at each of these points without rewarding them for the different rate of the growth of 
volume of the universe near each of the points [T6l [T71 HH [191 12] • 

Unfortunately, the results of the investigation using this set of the probability measures 
depend on initial conditions. More importantly, most of the physical entities which could be 
associated with "points" did not even exist before the beginning of inflation: protons did not 
exist, photons did not exist, galaxies did not exist. They appeared only after inflation, and 
their total number is proportional to the growth of volume during inflation. 

This gave rise to a set of different volume- weighted and/or pocket-weighted probability 
measures [TZl [TDl 1201 12U 1221 I2SI I2SI H31 123 I2SI 12^- The main problem with this 
approach is the embarrassment of riches: The total volume of the universe occupied by any 
particular vacuum state, integrated over the full history of the eternally inflating universe, is 
infinitely large. Thus we need to compare infinities, which is a very ambiguous task, with the 
answer depending on the choice of the cut-off procedure. 

Of course, this may simply mean that if life is at all possible in any of these vacua, it is 100% 
guaranteed to appear in any of such vacua. Therefore, instead of calculating the probability to 
be born in any of these vacua, one may try to find out in which family of the vacua do we live, 
according to our observations. Then we can take these data as an initial input for all subsequent 
calculations, and study conditional probabilities for the quantities which we have not measured 
yet [13J. This is a standard approach used by experimentalists who re-evaluate the probability 
of various outcomes of their future experiments on the basis of the latest experimental data. 

This approach is quite informative despite its limited nature. However, it is very tempting to 
go beyond this method and try to find a preferable way to introduce the cut-off and evaluate the 
relative probabilities to be born in different vacua. The simplest probability measure proposed 
in [HI EH [M [2D] introduces the cut-off at the upper bound of integration with respect to 
time. Unlike many other measures, this measure does not require different rules for different 
situations, it does not suffer from the so-called Boltzmann brain problem, and it is well suited 
for solving the cosmological constant problem [13J. Since this cut-off procedure depends on time 
in an exponentially expanding universe, it gives the results that can be exponentially sensitive 
to the choice of time parametrization JTQ1 20J. 

This observation stimulated the search for other probability measures which are not sensitive 
to the choice of time parametrization [221 [231 [2U [251 [26], [27]. However, the results obtained 
by using all of these measures remain sensitive to the choice of the cut-off. As the authors of 
this set of measures emphasized, their proposals allow a certain degree of flexibility which can 
be used to further specify their predictions; see e.g. [221 [21] ■ Therefore the search for the best 
probability measure still continues. 

In this paper we will take a fresh look at the simplest volume-weighted probability measure 
[HI [UJ [TUl [20] . We will identify the physical reason for the exponential sensitivity of the results 
to the choice of time parametrization [TOl i20j. After that, we will propose an improved version 
of the volume-weighted measure, which does not suffer from this problem, at least in the simple 
models that we have been able to analyze so far. 
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2 Volume-weighted distribution and bubble formation 



The methods that we are going to use in our investigation were described in our previous paper 
|13j ; here we will only briefly remind the reader of the necessary ingredients, and mention 
certain limits of applicability of these methods. 

In order to investigate a possible dependence of our results on the choice of the time slicing, 
we will will study slices of equal time with time measured in units H^~ X M^~^ . The choice 
(3 = 1 corresponds to standard comoving time measured by number of oscillations. The choice 
(3 = would mean that the time is measured in units of H~ l , or, equivalently, in the number 
of e-foldings. In such a case all parts of the universe expand at the same rate; this corresponds 
to the pseudo-comoving probability distribution of Ref. [13J. We are unaware of any physical 
interpretation of other values of (3\ in this paper we will consider (3 in the interval < (3 < 1 
and keep track of the dependence of our results on (3 in order to check whether the final results 
depend on time parametrization. 

Consider first the theory with a potential shown in Fig. HJ 
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Figure 1: A potential with three dS minima and two sinks. 

All dSi states are unstable and decay to other states dSj with the decay rates T^. Strictly 
speaking, one should multiply by an additional factor 4^ iff -4 [21]. However, in this paper 
we will assume that all factors are exponentially small and ignore all sub-exponential factors. 
As we will see in the next section, the main problems that we encounter appear already in this 
approximation. 

Let us consider, for definiteness, the part of the universe of size H^ 3 in de Sitter state 
corresponding to the upper local minimum of the potential in Fig. [TJ (The results of 
our approach do not depend on this particular choice.) All points inside the box initially are 
causally connected. Each point of dS2 will transit either to dSi or to dS$, from where it may 
either return back or go to a collapsing universe corresponding to AdS s or AdS z . In this paper 
we will assume that collapsing parts disappear into a singularity and can be ignored for all 
practical considerations. Therefore we will call AdS s and AdS z terminal vacua, or sinks. 
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Following [13] , one may begin with the investigation of this process in comoving coordinates, 
i.e. ignoring the expansion of the universe. To get a visual understanding of the process of 
bubble formation in comoving coordinates, one may paint red all of the parts of the original 
box corresponding the state (IS2, paint blue the parts in the state dS\, paint green the parts in 
the state (IS3, and paint black all points in any of the two sinks. In comoving coordinates, the 
size and the shape of the original box does not change during the cosmological evolution; all 
points remain at the same positions (same coordinates), but the values of the fields at each of 
these points (their colors) may change. 

In the absence of the sinks, the whole universe would become divided into red, blue and 
green parts, with the relative comoving volume of each of the parts proportional to e Si , where 
Sj is the entropy of each of the vacua. On the other hand, in the presence of the sinks, each 
particular points will eventually end up in the black. However, this does not mean that the 
whole universe is going to die. Indeed, under our assumption that Tji are exponentially small, 
the total volume remaining in each of the dS vacua will continue growing exponentially. This is 
a crucial feature of the eternal inflation scenario, but it is hard to see it in comoving coordinates, 
ignoring the growth of volume. 

Here we have two different possibilities. If we associate ourselves with some particular point 
which existed at the moment of creation of the universe, then we can study the most probable 
evolution of this point during the time until it dies in the sink. One could hope that by doing 
so we will find the most probable trajectory which could bring us to the present state. This 
would correspond to the evaluation of comoving coordinate probabilities [T6l [T71 [T8| [T9]. 

However, this proposal is rather problematic. Let us assume, for example, that when our 
universe was born, it was a compact universe of a Planckian size. According to [211 E2] , quantum 
creation of such universes seems most probable. This result was recently used in one of the 
papers advocating the usefulness of comoving probabilities [30J. But Planckian-size universes 
in the beginning could contain just a few physically distinguishable points: We do not know 
any way to distinguish between points separated by less than the Planckian distance. 

Thus we cannot associate ourselves with anything existing at the beginning of inflation: The 
total number of physically distinguishable points in the universe at the moment of its creation in 
the model discussed above is much smaller that the number of people living now on the Earth. 
A similar conclusion is true for most of the inflationary models independently of their initial 
state: If we consider any two observers living now on the Earth, the distance between them 
at the beginning of inflation was smaller than Planckian. Moreover, all elementary particles in 
our bodies did not exist before inflation, and even the galaxy where we live did not exist before 
inflation. Therefore it is hard to imagine that the behavior of the physical conditions at a given 
initial point, ignoring production of many new points during the cosmological evolution, could 
have any anthropic significance. In a certain sense, the eternal existence of an infinite number 
of points where life may or may not emerge is just a bad mathematical idealization, which is 
especially misleading when applied to inflationary cosmology. 

One of the main roles of inflation is to create exponentially large number of elementary 
particles in each part of an inflationary universe which originally occupied an extremely small 
causally connected region of the size Oi^H^ 1 ). But in order to study the new parts of the 
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universe eternally created by inflation one should pay attention to the exponential growth of 
the volume of the universe. One of the ways to do it is to study the volume-weighted probability 
distribution introduced in [HJ [TQl [20]. In application to the scenario with transitions between 
the metastable dS vacua, the set of equations for the probability density Pj for the potential 
with three different dS minima and two AdS minima was given in 



p 1 = -p 1 (r ls + r 12 ) + p 2 r 21 + 3H?p 1 , (1) 

p 2 = -p 2 (r 2s + r 21 + r 23 ) + p 1 r 12 + p 3 r 32 + 3p^p 2 , (2) 
p 3 = -p 3 (r 3z + r 32 ) + p 2 r 23 + 3h$p 3 . (3) 

The first few terms in these equations describe the transitions between different vacua. The last 
terms describe the growth of volume of de Sitter spaces dSi. These equations intuitively are 
very simple, but one should keep in mind that their derivation involved several approximations 

m 

We will be interested in the case where is much greater than all other parameters in 
these equations. In this case P 2 obeys a simple equation 

P 2 = 3H$P 2 , (4) 

i.e. in the first approximation P 2 does not depend on Pi, P 3 : 

P 2 = P 2 (0)e 3 ^* . (5) 

Because of the fast growth of P 2 , the terms P 2 r 2i eventually become the leading terms in the 
equations for Pj, for all i ^ 2: 

P = p 2 r 2i . (6) 

The solution of this equation is 

m) = ^^ (l-e-^') . (7) 

Within a very short time O^H^ 13 ) the solution approaches the asymptotic regime 

P 2 (t)r 2l _p 2 (o)r 2j w, t 



The ratio of Pj is given by 



3P 2 P 3P 2 P 

where i,j ^ 2. Note that this result does not depend on the time parametrization described by 
the factor (3 in the approximation specified above. 

Applying these results to anthropic considerations requires several additional steps. The 
functions Pj describe the total volume of all parts of the universe in a given vacuum state, 
produced from the initial domain. However, when the bubbles of a new phase expand, their 
interior eventually becomes an empty dS space. If we are usual observers born after reheating 
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of the inflationary universe, then one may argue that the probability to be born in the bubble 
dSi during any given time interval is proportional not to the total volume of de Sitter bubbles, 
which is described by Pj, but to its part produced within some finite time after the bubble 
formation, e.g. in the interval from 1 to 100 billion years. In other words, the probability of life 
is proportional not to the total volume, but to the new volume, which switches our attention 
form the total probability Pj to the incoming probability charge Qi and the incoming probability 
current Qi [20, HSJ Q35]. When we say 'incoming,' we emphasize that we are interested only in 
the growth of the volume of the new parts of the bubbles, but not in the speed of their decay. 
This growth occurs due to the new bubble production, and due to the expansion of the walls 
of the existing bubbles. In the case where is much greater than all other Pf , the leading 
part of this process occurs due to the new bubble production, and equations for Qi have the 
following simple form, coinciding with the form of the leading terms in the equations for Pf 



In this case, the ratio of the probability currents from the upper minimum to the lower 
minima will be given by 



for i,j ytz 2. As before, this result does not depend on the time parametrization described by 
the factor (3 in the approximation specified above. 

Thus, the predictions obtained using the volume-weighted probability distribution do not 
depend on the choice of time slicing for a wide variety of time parametrizations. This depen- 
dence appears only for exponentially small f3, i.e. in the limit when one completely ignores 
the expansion of the universe. However, in this limit our equations should be written more 
precisely. 

Our results do not depend on initial conditions either. Indeed, if we assume that initially 
the field was positioned not in the highest minimum, the universe will expand at a slower rate 
until the jump to the highest minimum occurs. Starting from this time, the volume of the part 
of the universe in the highest minimum (and the total volume of its descendants produced by its 
decay) will grow much faster than all other parts. Eventually this will result in the stationary 
regime discussed above, with o i? - = = 3± = S. — Im. 

■> Uij Qj fj 1 2j 

Let us describe our approximation again. First of all, we assumed that is much greater 
than all decay rates T^. This is a legitimate approximation as long as we ignore all sub- 
exponential factors. More importantly, we assumed that is much greater than all other 
Pf . If P 2 > Hi for i ^ 2, the condition H$ > Pf is violated only for (3 < log" 1 g < 1. 
We assumed also that as soon as bubbles of the vacuum dSi form, their interiors start growing 
exponentially, as e 3H * *, and starting from the moment of their formation, decays to other vacua 
with the decay rate T^- become possible. 



Qi = P 2 r 



21 



(10) 




(11) 

(12) 




(13) 
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This approximation seems reasonable, but one should remember that it is just an approxi- 
mation. For example, if one restores the sub-exponential term in the definition of T^, 
one will encounter the sub-exponential dependence of the results on (3. It may happen that this 
dependence will disappear when one solves all other parts of the problem without making any 
approximations. We will not discuss this issue here since we have a much more complicated 
problem to solve: Once we add to this picture a stage of slow-roll inflation, the dependence on 
f3 becomes exponential [20J. In the next section we will discuss this problem, identify its source, 
and propose its possible resolution. 

3 Volume-weighted measure for slow-roll inflation 

Let us now consider inflation in a theory with the potential shown in Fig. [2j which describes 
both the tunneling and the slow-roll. Inflation begins while the field is at the top of the 
potential, = 3 . Then the field tunnels either to 04 or to 02, and experiences a slow- roll 
regime until it reaches either 5 or 1; which corresponds to the end of inflation boundary. 
Once this happens, the field instantly falls into the minimum nearby, and the corresponding 
part of the universe becomes hot. 



V 




Figure 2: Tunneling from the false vacuum to two different slow roll inflation regimes. 

The first part of this process can be described by the equations similar to the ones studied in 
the previous section, with an obvious result: At H^t ^> 1 the probability distributions rapidly 
approach an asymptotic regime 

P 3 (t) ~ P 3 (0) e 3H "\ Pa = Qa = P 3 (0) r 34 e 3H 3\ P 2 = Q 2 = P 3 (0) r 32 e 3 ^ . (14) 

Just as before, the ratio Q4/Q2 = Q4/Q2 = P4/P2 = ^4/1^32 does not depend on time t and 
on the time parametrization (3, which is encouraging. 
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However, let us now discuss the subsequent slow-roll stage. This stage cannot be described 
by the simple equations used in the previous section. Equations for the tunneling allowed 
almost instant response for the exponential growth of the probability distribution P 2 : Once the 
total volume of the universe at the top of the potential starts growing exponentially as e 3H2t , 
this growth almost instantly results in the proportional growth of all other Pj. The situation 
during the slow- roll regime is qualitatively different: The field cannot tunnel anywhere, and 
the universe cannot reheat until the slow-roll regime is over. One can describe this effect using 
the volume- weighed Fokker-Planck equations [10], or using a simplified slow-roll description 
ignoring quantum diffusion of the scalar field on its way down. Here we will use this simplified 
approach; it will be quite adequate for our purposes. 

After tunneling to the point 04, the field slowly rolls to the left during some time £45 until 
it reaches the point 05, where inflation ends. Suppose for simplicity that during this time the 
universe inflates with the Hubble constant H 45 « const. Similarly, the field rolls to the right 
from 02 during some time t 2 i until it reaches the point 0i where inflation ends. During this time 
the universe inflates with the Hubble constant H21 ~ const. We will assume, for definiteness, 
that £21 > £45- 

Because of the rapid growth of the upper dS space, the total volume of all parts of the 
universe eventually approaches a stationary regime, when the total volume of all parts of the 
universe will grow approximately as e 3H ' 2t . This is a very counter-intuitive result, which we 
already discussed in the previous section, but now we need to understand it in a much more 
detailed way. 

Let us consider this process step by step. Suppose the tunneling begins at the moment 
t = 0. The number of bubbles with the field 02 at their center will be the same as the number 
of the bubbles with the field 04 at their center. (Note that in the coordinate system that we are 
using the bubbles have finite size, and the initial value of the field at their center is determined 
by the CDL tunneling trajectory.) But during the time interval £21 after the tunneling there 
will be no bubbles with the field 0i at their center. When they appear, their total volume 
will be equal to the total volume of the first bubbles with the field 2 , multiplied by e 3N21 , 
where N 2 i = H 21 t 2 i is the number of e-folds of the slow-roll inflation on the right branch of 
the potential. Note that this factor does not depend on (3. The factor e 3N21 describes the total 
growth of volume during the slow-roll inflation; it is gauge invariant. 

Starting from t = £21, the total volume occupied by the field 0i grows as e 3H3 ^~ t21 ^ : 
P(0i,t) = for t < t 21 , 

p(0 1; t) = p 3 ( ) r 32 e 37V21 e 3if 3 for t > t 21 . (15) 

Similarly, for the left branch we have 

P(05,t) = fort <U 5 , 

P(0 5 ,t) = P 3 (0) r 34 e 3Ni5 e 3H * ( '-' 45) for t > t 45 . (16) 

If we literally follow the rules of (TD|, we should divide P(0i, t) by P(0s, t) at large t, to obtain 
the following result: 

P(0 1; t) _ r 32 e 3 ^ HlM 

P(05,*) r 34 e3^ 5 ■ ^ > 
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Here 

At = t 45 - t 21 . (18) 

This result is exponentially sensitive to the time parametrization factor (3 [10J. This was the 
main fact that forced many authors to look for a better probability measure for eternal inflation. 

However, in the simple example considered above one can easily identify the source of the 
problem. We were able to find a time-independent result for the ratio of the two probabilities 
because both of the expressions P(4>i, t) and P(4>5, t) asymptotically approach the same station- 
ary regime; they grow at the same (exponential) rate e 3 ^ 2 *. But they approach the stationary 
regime at two different moments of time because of the time delay At = t 45 — t 21 related to 
the slow roll regime. 
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Figure 3: Large green bubbles correspond to the bubbles containing the field Small red 
bubbles contain field 05. The chain reaction of the bubble production starts at the time t\ 2 for 
the 'tree' of the bubbles of the field 0i, and at the time £45 for the 'tree' of the bubbles of the 
field 05. After this process starts, the number of the bubbles belonging to each 'tree' grows at 
the same rate, as e 3H2T , where r is the time starting from the moment when the first bubble 
of each type appears. Instead of comparing the number of the bubbles at the same moment 
of time t [TO] , which does not make much sense for t < max{t 2 i, t i5 }, we propose to reset the 
time for each type of the bubbles to the moment r = 0, when their number starts growing at 
the same rate, and compare them at the same time r. In other words, one should compare 
apples to apples, instead of comparing apples to the trunks of the trees. The results of this 
comparison in our simple model do not depend on initial conditions, on the time r, or on the 
choice of time parametrization. 

In the context of this scenario, it does not make much sense to compare two time-dependent 
quantities before they have reached the asymptotic stationary regimes. This problem was 
hidden in the original proposal of [101120] because we studied there the limit t — > oo, where this 
issue could not be clearly identified. Therefore we believe that our original proposal to compare 
different processes at the same moment of time, starting from the beginning of inflation, is 
incomplete and requires a certain modification. Indeed, the probability distribution proposed 
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in [TO], EH] was very far from stationary at the beginning of the process, when £ < maxfazi, £45}, 
see Fig. [31 The simplest modification is to compare the processes in different vacua (or at 
different values of the scalar fields, or temperature, or density) starting at the time when the 
stationary regime becomes established in these states. The corresponding moments of time will 
be different for the vacua of each particular type. We will call the improved probability measure 
stationary, to distinguish it from the asymptotically stationary measure proposed in [TO"! [20]. 

This means that we should reset our clock, ignoring the time delay related to the different 
onset of stationarity for the parts of the universe with different properties. In the example 
considered above, this means that we should compare the volume at the point <pi at the time 
r = t — £ 2 i with the volume at the point </> 5 at the time r = £ — £ 45 . This leads to the improved 
expression for the ratio of the two probabilities: 

Willi = ^ ^ m 

This is the main result of our paper. This result is gauge-invariant, it does not depend on (3, 
unless one takes /3 < 1, which violates the validity of our approximation. This result has a very 
simple intuitive meaning: It shows that the fraction of the volume of the universe in a given 
state is proportional to the probability T 3i to go in a given direction, and to the total growth 
of volume along the corresponding inflationary trajectory. 

Note, that in the derivation of this result we used an assumption that the evolution of the 
probability distribution is dominated by the expansion of the domain with the highest energy 
density. This is indeed the case for all but exponentially small (3. However, this condition is 
not satisfied in the limit (3 = 0, which would correspond to the comoving or pseudo-comoving 
probability distribution [13] . In other words, the results of our calculations in the simple model 
discussed above do not depend on time parametrization for all (3 except for the limiting case of 
(3 = 0. 

This resembles the discussion of the gauge-noninvariance of the calculation of the critical 
temperature in the unitary gauge as compared to (almost) all other gauges [33j. The difference 
appears there because of the use of non-invariant quantities, such as the scalar potential, and 
also because of the bad infrared behavior of the higher-loop corrections near the critical point 
in the unitary gauge. Eventually we have learned how to calculate the critical temperature in a 
gauge-invariant way, but the corresponding procedure is quite involved, so the original simple 
but slightly non-invariant methods still remain valuable and widely used. 

A possible reason for the particular behavior of the probabilities calculated in the limit (3 — > 
is that for all sufficiently large (3 > 0, the time that is required for reaching the stationarity 
regime depends on the parameters of the theory logarithmically, the qualitative picture of the 
flow from the upper minimum remains correct, and the calculation of the delay A£ does not 
change. Meanwhile for j3 < the probability current would flow from the lowest dS vacuum, 
and the delay time A£ would be exponentially large. Using the analogy with the theory of phase 
transitions, one would argue that one should calculate the probabilities by our method only if 
13 is not too small; the calculations in the limit j3 = using our method become unreliable. 
Admittedly, this is just an argument; this issue requires a more detailed analysis. 

Let us formulate our results in a different form. Suppose we are interested in the relative 
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volumes of the parts of the universe with two different temperatures T\ and T 2 in the two 
different vacua. The generalization of our previous result is obvious: 

HTut) r 32 V(TQ 

P(T 2 ,r) r 34 V(T 2 )- {U) 

Here V(T 2 ) is the total growth of volume of the universe from the moment of tunneling to the 
point 2 to the moment when inflation ends and the temperature in the corresponding part of 
the universe becomes equal to T\. Similarly, V(T 2 ) is the total growth of volume of the universe 
from the moment of tunneling to the point 04 to the moment when inflation ends and the 
temperature in the corresponding part of the universe becomes equal to T 2 . 

If we would like to compare the volumes for the same vacuum state, at different tempera- 
tures, we would get a trivial result 

PgVr) _V(T 1 ) 

P(T 2 ,r) V(T 2 ) ' 1 ) 



It is instructive to compare this result to the result which we would obtain using the old 
measure of PUJ [2U] prior to its improvement: 

P(Tur) _ V(T\) HlAt 

P(T 2 , r) " V(T 2 ) ' [ZZ) 

where At = t(T 2 ) —t(Ti), and t{T) is the time when the universe cools down to the temperature 
T. Clearly, our old measure would imply exponential dominance of the high temperature regions 
|34j, which was the essence of the so-called youngness paradox [35]. However, once we reset 
our clock so that the comparison begins starting from the time when the stationary regime is 
established, the dangerous factor e H * At leading to the youngness paradox disappears. 

One should note that this was indeed a paradox rather than a problem, since the volume 
distribution obtained in [10] was not supposed to be used in anthropic considerations. Mean- 
while the volume-weighted anthropic probability measure, which was first proposed in [20], uses 
the methods developed in |l0j . but calculates the probability currents through the hypersurface 
of a given temperature/density/scalar field. With the probability measure proposed in [20J, the 
youngness paradox does not appear [13]. (This fact remained unnoticed in many subsequent 
discussions of the volume-weighted probability measure, see e.g. [3UE5].) However, it is quite 
encouraging that with our improved prescription for volume-weighted probabilities, the young- 
ness paradox disappears altogether, even if one decides to compare completely different regions 
of the universe, without restricting this comparison to the hypersurface of the end of inflation, 
or of a given temperature or density. 
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4 Discussion 



In this paper we proposed an improved volume-weighted probability distribution for the mul- 
tiverse. We performed our investigation ignoring sub-leading factors. For the simple model 
described in Sections |2] and El our approach produced gauge- invariant results, thus avoiding the 
problem plaguing the probability measure proposed in pH [20] • The results obtained by using 
our improved volume-weighted measure, which we called stationary, have a simple intuitive 
interpretation, which shows that probability is proportional to the growth of volume during the 
slow-roll inflation. 

The main reason for the gauge invariance of our result can be understood as follows. In any 
time parametrization, the probability distribution becomes time-independent once the asymp- 
totic stationary regime is reached. In our model the post-inflationary stationarity is established 
when the field inside the bubbles produced by the tunneling from the state dS$ reaches the 
end of inflation (points 1 and 5). Since the definition of the end of inflation is gauge-invariant, 
the probability distribution at the corresponding hypersurface is also gauge- invariant, and it 
remains gauge invariant at all subsequent time slices because of the stationarity. 

It is interesting to compare the results of our approach with the results which could be ob- 
tained by other methods. At the first glance, our approach may resemble the method developed 
in [22] , where the cut-off procedure depended on the fraction of the thermalized comoving vol- 
ume. However, these methods are quite different. In particular, our results, unlike the results 
of [22], do not depend depend on initial conditions. 

Our main result for the probability distribution ffT9l in the simple model with the potential 
shown in Fig. [2] coincides with the result obtained by Vilenkin et al using a different probability 
measure developed in [24, 25J. However, this coincidence disappears when one studies more 
complicated models. For example, one may consider a potential with an additional de Sitter 
minimum dS^, somewhere on the way from the point 2 to point 1. As one may anticipate from 
the results obtained in the previous section, this should lead to an additional suppression of 
the probability to live in the minimum near the point 1 by the exponentially small factor r 61 . 
The meaning of this result is very simple: The exponentially growing probability current flows 
from the upper vacuum in the direction of the least resistance. When we add an extra barrier, 
the corresponding current slows down. This conclusion can be easily verified by solving the 
corresponding system of equations. 

In the similar situation studied in [25] the authors identified two possible results of the 
introduction of the additional minimum. The first of the two probability measures discussed 
in [25] predicts that the insertion of an additional metastable minimum dS§ does not affect 
the ratio of the probabilities to live in the right minimum as compared to the left minimum. 
Meanwhile the second of these two measures, which was advocated in [25J, suggests that the 
insertion of the additional minimum makes the relative probability to live in the left minimum 
vanishingly small. This is directly opposite to the effect expected in our approach. Note that 
this difference appears even in the models without the slow-roll regime. 

This conclusion, however, requires a more detailed investigation. The calculation of the 
time delay for the simple model shown in Fig. [2] was quite unambiguous. In more complicated 
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models with many minima it may be much more difficult to identify the time when the stationary 
regime is established. A small error in determining this time may translate into a very large 
error in the relation between the probabilities, proportional to e 3 ^ 2 At . We did not find any such 
terms in the models without the slow-roll regime, but it is quite possible that a more accurate 
investigation will reveal them. 

The situation may become even more interesting and complicated if one considers Hawking- 
Moss tunneling, or includes the stages of eternal slow-roll inflation. This will require an inves- 
tigation of the branching diffusion equations along the lines of [TUl [20] . Thus the calculation of 
the probabilities in the landscape can be a rather challenging task. But nobody expected an in- 
vestigation of the eternally growing inflationary multiverse to be simple. Instead of complaining 
about the hardship of life, let us try to concentrate on the positive news. 

There are many different proposals for the calculation of the volume- weighted/pocket- 
weighted probabilities during eternal inflation. Now we reached the point where the most 
popular proposals for the volume-weighted or pocket-weighted probability measure give iden- 
tical results at least for the simple model studied in our paper. This result is quite simple 
and intuitively clear: The prior probability to live in a given dS state is proportional to the 
probability of its formation from the dS vacuum with the highest energy density, multiplied by 
the growth of volume during the slow-roll inflation. 
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